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Functionality in single-molecule devices: Model calculations and applications 
of the inelastic electron tunneling signal in molecular junctions 
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We analyze how functionality could be obtained within single-molecule devices by using a combination of non- 
equilibrium Green's functions and ab-initio calculations to study the inelastic transport properties of single- 
molecule junctions. First we apply a full non-equilibrium Green's function technique to a model system with 
electron-vibration coupling. We show that the features in the inelastic electron tunneling spectra (IETS) of the 
molecular junctions are virtually independent of the nature of the molecule-lead contacts. Since the contacts 
are not easily reproducible from one device to another, this is a very useful property. The IETS signal is 
much more robust versus modifications at the contacts and hence can be used to build functional nanodevices. 
Second, we consider a realistic model of a organic conjugated molecule. We use ab initio calculations to study 
how the vibronic properties of the molecule can be controlled by an external electric field which acts as a 
gate voltage. The control, through the gate voltage, of the vibron frequencies and (more importantly) of the 
electron-vibron coupling enables the construction of functionality: non-linear amplification and/or switching 
is obtained from the IETS signal within a single-molecule device. 



I. INTRODUCTION 

Single-molecule electronics has shown significant 
progress in recent years. Aviram and Ratner first pro- 
posed in the 1970s electronic devices where a single 
molecule functions as the active element^. Since then 
a variety of interesting effects has been observed in the 
transport properties of single molecules, including recti- 
fication, negative differential resistance, and switching. 
The molecules usually employed are organic conjugated 
molecules with strongly delocalized 7r-electrons along the 
molecular backbones (and eventually functional chemi- 
cal groups acting as electron donors or acceptors), and 
are natural candidates for bottom-up assembly of elec- 
tronic devices. Chemical synthesis is a massively parallel 
method to create identical molecules by the mole. Their 
structures and electronic properties can be tailored vir- 
tually at will and can be characterized by several spec- 
troscopic and structural probe techniques. Thus they are 
ideal building blocks for electronic devices. 

However, to make even the simplest molecular elec- 
tronic device — a single molecule connected between two 
electrodes — presents several practical difficulties. The 
fabrication of the nanojunction is a challenge because 
of the difference in scales between the molecule, the elec- 
trodes, and the small gap between them. Moreover, it is 
even more difficult to verify that the intended molecule 
is indeed in the junction, that it is well connected, and 
that it is oriented in the expected manner with respect 
to the electrodes. Several methods are used to create 
such nanojunctions (e.g. mechanical break junctions 2 ^, 
crossed wires^, and scanning probe microscopy^). Each 
technique has its own advantages. For example, some of 
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them can be cooled to low temperature, others provide a 
good control of the electrode spacing, or provide images 
of the contact area. However, the most crucial and as yet 
uncontrolled point is the nature of the contacts between 
the molecule and the electrodes. This is pivotal, since the 
transport properties of nanojunctions do not depend just 
on the intrinsic properties of the molecules, but also on 
the nature of the contacts. Although there is an increas- 
ing amount of research into the control of these contacts 
by tailoring the end groups of the molecule anchoring on 
the surface^— , there is still a great deal of uncertainty 
about the specific nature of the contacts. These change 
from device to device, or even fluctuate during a single 
experiment, because of thermal effects (atom migration 
at the electrode surface) or of current-induced effects. 
Hence it is extremely difficult to produce molecule junc- 
tions with reproducible properties. 

Other properties of current-carrying molecular junc- 
tions are also important to understand, such as the na- 
ture and the effects of electron-vibron coupling on the 
molecule. Such a coupling can modify the transport 
properties of the junctions at a particular threshold volt- 
age and is also responsible for heating effects in the junc- 
tions. Several groups have already performed inelas- 
tic electron tunneling spectroscopy (IETS) on molecular 
junctions to obtain a signature from the vibration modes 
of the molecules in the nanojunction^r— . 

The realization of a true single-molecule device, with 
three terminals, has been achieved recently^. The trans- 
port properties through the source and drain are directly 
modulated by an external gate voltage as in a conven- 
tional field-effect transistor. However, in the molecular 
transistor, the gate voltage modifies the molecular orbital 
energies. In Ref. [ltl the authors also studied the effects 
of the gate voltage and temperature on the IETS signal. 
They found that depending on the nature of the molecule, 
the IETS signal is nearly independent of the gate volt- 
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age (for a molecule with a cr-saturated alkyl backbone), 
however significant modifications of the IETS features' 
intensity and lineshape exist for 7r-conjugated molecules 
with an aromatic ring. Their results indicate that both 
conductance and IETS signal can be modified and con- 
trolled by a gate voltage. 

Full ab-initio calculations for a realistic single-molecule 
with three terminals (source, drain and gate) are ex- 
tremely computationally demanding, especially for self- 
consistent calculations and/or calculations for realistic 
non-equilibrium conditions. Non-equilibrium transport 
properties have been calculated for realistic systems of 
single-molecules connected to two terminals by using 
density-functional theory22r— . Calculations in the pres- 
ence of a third gate electrode are rare^ir— and usually 
use model Hamiltonians (i.e. tight-binding or extended 
Hiickel). In this paper we address the problem of func- 
tionality in three-terminal single-molecule devices by us- 
ing a combined (two step) theoretical framework. We 
discuss the use of both the electronic and the vibronic 
properties of the molecule to achieve a device with func- 
tionality. 

The first step consists of a detailed analysis of the full 
non-equilibrium transport properties of a model system. 
We study the effects of the contacts between the molecule 
and the source and the drain on the conductance and the 
IETS signal. The presence of the gate voltage is incorpo- 
rated in the position of the molecular levels and generates 
transitions between different transport regimes (resonant 
and off- resonant). We use a full non-equilibrium Green's 
function (NEGF) technique^i 3 ^ 3 -^— . In order to explore 
a wide range of experimental contact geometries, we first 
study a model system, in which two parameters (the con- 
tact strength and the corresponding fractions of potential 
drop) 40 ' 59 characterize the molecule-lead coupling. We 
show that the IETS signal is much less sensitive to the 
nature of the contact than the conductance itself. 

The second step consists of applying the principle 
found in the first step to a more realistic model. We 
study the influence of an external electric field on the 
electronic and vibronic properties of a realistic molecular 
system based an ethynylphenyl-based backbone by using 
ab-initio calculations^ -. Having shown, for the model 
system, that the IETS signal is nearly independent of 
the nature of contacts, we concentrate on an isolated 
molecule and extract the relevant physical quantities (vi- 
bration frequencies and electron-vibron coupling matrix 
elements) in the presence of the electric field which sim- 
ulates the presence of a gate electrode. 

Finally we discuss and show the possibility of using 
the IETS signal, rather than the current or the con- 
ductance, to form functional nanodevices with switching, 
non-linear amplification or sensor functionality by using 
our ab-initio calculations for the conjugated molecule. 



II. NON-EQUILIBRIUM TRANSPORT AND THE 
EFFECT OF THE CONTACTS 

We use a NEGF technique to calculate self-consistently 
the full non-equilibrium inelastic properties of a molecu- 
lar junction. Using a model system to reduce these cal- 
culations to a tractable size, we concentrate on a sin- 
gle molecular level coupled to a single vibrational mode. 
A full description of our methodology is provided in 
Refs. and 113 ■ 

This approach allows us to study how variations at the 
contacts modify the transport properties of the junction. 
These include experimentally uncontrollable modifica- 
tions of the molecule-lead coupling due to the experimen- 
tal environment — such as thermal fluctuation, diffusion 
of atoms at the surface of the leads, variation of the gap 
of the nanojunction, the presence of impurity molecules 
around the nanojunction, etc. The consequence of these 
modifications, and the reason they limit the reproducibil- 
ity of functionalized junctions, is to change the geometry 
of the contacts, and therefore the strength of electronic 
coupling between the molecule and the leads, and the 
corresponding potential drops at the contacts. 

In our model, the central region of the nanojunction, 
i.e. the molecule, is described by a simple electron-vibron 
coupling Hamiltonian 

He = Eod'd + huoa'a + 7o(a^ + a)dJd, (1) 

where one electronic molecular level £o and one vibra- 
tion mode with energy ujq are coupled together via the 
coupling constant 70. The electron and phonon cre- 
ation/annihilation operators are S jd and a^/a respec- 
tively. 

The central region is then coupled to two non- 
interacting Fermi seas each at their own equilibrium, 
characterized by the left and right Fermi levels /i£ and 
fin, via two hopping matrix elements toL and ton which 
represent the strength of electronic coupling between the 
molecule and the leads. 

All the properties of the nano junctions (electronic 
density, spectral functions, current density) are deter- 
mined from the knowledge of the NEGF of the central 
regio n 34 ' 35 ' ■ For example, the retarded Green's func- 
tion of the central region is given by 

G» = [u>-e - S7 L M + YT R (w) + K- V -Mr\ (2) 

where Y7 a is the self-energy arising from the a = L,R lead 
and Sg_ vib is the self-energy arising form the electron- 
vibron interaction. The latter is calculated from Feyn- 
man diagram expansion of the interactional. In this sec- 
tion, we calculate £ e _ v ;b from the lowest order expansion 
of the electron-vibron interaction (first Born or Hartree- 
Fock level of approximation) 34 i 35 . The calculations are 
performed self-consistently^ 4 ^. 

Within our NEGF model, the electrostatics are not 
solved self-consistently with the non-equilibrium elec- 
tron charge density, and we thus have another degree 
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FIG. 1. Schematic representation of our model system. Case 
(I) shows a symmetric coupling at the contacts, with tf )L — 
tg R . Cases (II) to (IV) are asymmetric junctions. In Case (II) 
the total gap is kept constant, but there is a variation in the 
position of the central region, i.e. t f + tg R = t§ L + t§ R . In 
Case (III) the total gap is compressed, tg£ + t^R < ^ol + ^aR) 
while in Case (IV) it is expanded, tg£ + tg R > tf )L + t^ R . 

of freedom for the potential drops at the contacts. At 
equilibrium, the whole system has a single and well- 
defined Fermi level /x oq . Out of equilibrium, a finite 
bias is applied across the junction, giving Fermi levels 
fJ>L,R = M cq + VL,RfV. Following Ref. |59j . the frac- 
tion of electrostatic potential drop at the left contact is 
r]L = +t]v and 77^ = — (1 — rjv) at the right contact, with 
i]L — T)r — eV and rjv <E [0, 1], as shown in Fig. [TJ 

The parameter rjv characterizing the potential drop 
at the contact is related to the other parameters of the 
junction, as we can see from the following: Let us con- 
sider three typical cases. For the first, when the coupling 
at the contacts is symmetric, i.e. when £ql = ioi?., it is 
reasonable to assume a symmetric potential drop, i.e. 
rjv = 1/2. In the second case, the coupling is very strong 
on one side, for example ioi ^ toR, the left Fermi level 
fiL is pinned to its equilibrium value ^ cq and there is a 
large tunneling barrier at the right contact where all the 
potential drop occurs, i.e. ji R — fi cq + eV. Finally, we 
have the opposite case, when £ol ^ tQRi the right Fermi 
level flu is pinned to /x eq all the potential drop occurs at 
the left contact, i.e. hl = A*° q + e V ■ 

For intermediate cases, it is therefore reasonable to as- 
sume a lowest-order linear dependence of the parameter 
rjv with the hopping integral parameters ioL.i?, and hence 



we choose the following relation rjv — tou/(toL + 
We now have a relation between the strength of elec- 
tronic coupling between the molecule and the leads and 
the corresponding fractions of potential drop. In reality 
the hoping integrals depend exponentially on the distance 
between atoms and the dependence of the fractions of po- 
tential drop on such a distance obeys a different relation. 
In principle, the later depends on the dielectric properties 
of the fully connected molecular junctions. Example of 
fractions of potential drop rj calculated for realistic molec- 
ular junctions, with symmetric and asymmetric coupling, 
are given in Ref |6l| . 

In the following, we will consider four different cases, 
as depicted in Fig. [TJ corresponding to different possible 
modifications of a symmetric nanojunction, and calcu- 
late their transport properties. However, in all the four 
cases, the fractions of potential drop will be close to the 
symmetric case. 

We concentrate on the regime corresponding to tunnel- 
ing through the HOMO-LUMO gap of an organic molecular 
junction, where the molecular level is well above (or be- 
low) the equilibrium Fermi level /i cq ). This is a typical 
behavior for a semiconducting-like molecule sandwiched 
between two electrodes when the gate voltage is small. 

The calculations for the NEGF model system can be 
reduced to unitless parameters, i.e. normalized by the 
hopping integral for the semi-infinite one-dimensional 
leads. We present the most relevant results of our study 
for the following set of parameters: eo = 1.5, wo = 
0.4, 70 = 0.35 and two sets of coupling to the leads: 
medium coupling with t$ L = tQ R = 0.27 and weak cou- 
pling with £g L r — 0.15. A wide range of parameters 
has been explored, and this choice is the closest to those 
calculated from first principles (see below). 

The fluctuations in the nanojunction introduce vari- 
ations of the hopping matrix elements shown in Fig. [1] 
from a symmetric coupling junction, we obtain an asym- 
metric junction when the position of the molecule inside 
the gap is varied by a small amount (10%). Hence the 
left and right hopping integrals are tfif = t^ L + 10% and 
*or = toR ~ 10% respectively — case (II). We also con- 
sider two cases for which the gap is modified according 
to t^ R — t^f for a compression of the gap — case (III) — 
and as t^ R — t$ R for a expansion of the gap — case (IV). 
The potential drop parameter rjv that determines the 
fractions of potential drop at each contact is calculated 
from rjv = toi?7(ioL + toR) as we have explained above. 
For our set of parameters, we get rjv — 0.50 for case (I), 
■qv = 0.45 for case (II), rjv = 0.5238 for case (III) and 
7j V = 0.4737 for case (IV). 

The corresponding conductance curves are shown in 
Fig. [21 and show one main conductance peak for each 
configuration of the junction as expected. The con- 
ductance peak corresponds to a resonant transmission 
through the main electronic level of the central region, 
renormalized by the electron-vibron coupling, i.e. a peak 
at s» eo ~ eo — lo/^o- It is clear from Fig. [2] that the bias 
position of the peak depends strongly on the value of the 
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potential drop factor r\y when all the other parameters of 
the central region are kept unchanged. The width of the 
peak, proportional to t^ L + tQ R , is not greatly affected by 
the 10% fluctuation of the contacts for the values of the 
parameters we have chosen. 




0.8 1 1.2 
V[V] 



FIG. 2. Dynamical conductance G(V) = dl/dV versus ap- 
plied bias for the four cases depicted in fig. [1] symmetric 
(I) and asymmetric cases (II, III, IV). The top panel shows 
medium coupling to the leads (tg L — 0.27) and the bottom 
panel weak coupling (tg L — 0.17). Both sets of calculations 
show that, even for small variations of the potential drop pa- 
rameter rjv, the position of the main conductance peak is 
strongly dependent on the fraction of potential drop at the 
contacts. 



Results for the IETS signal for the same parameters 
are shown in Fig. [3j where we clearly observe a rather 
different behavior in the conductance peaks. Here we see 
two separate features — the first is a peak at the vibra- 
tion energy ojq, as we would expect for the off-resonant 
regim o 35 ' i 40 . The amplitude of this feature is propor- 
tional to 7q while the width is dependent on the other 
parameters of the junction o 56 i 62 . At higher biases we ob- 
serve the conductance peak, which occurs at the polaron- 
shifted electronic level £q. Increasing the coupling to 
the leads yields an effective decrease of the conductance 
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FIG. 3. IETS signal d 2 I/dV 2 normalized by the conductance 
G(V) for the four cases depicted in fig. [1] The top panel 
shows data for weak coupling to the leads (tg L = 0.17), the 
middle and bottom panels for medium coupling (tg L = 0.17). 
Both sets of calculations show that the position of the IETS 
feature at V — luq = 0.4 is not dependent on the fraction 
of potential drop at the contacts. The bottom panels show 
different normalizations of the IETS as well as the bare IETS 
signal, as different normalization conventions are often used in 
the experiments. There is no need to show (d 2 1 / dV 2 ) / (I /V) 
since in this bias regime the current is linear G ~ I/V. 



peaks accompanied by an increase in the width, as ex- 
pected. This leads to better contrast in IETS between 
the pure inelastic feature at uiq and features associated 
with elastic and inelastic resonant tunneling: the ampli- 
tude of the inelastic feature is much less dependent on 
the coupling to the leads. 
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It is also clear from Fig. [3] that the position in bias of 
the IETS features around V — uj does not depend at 
all on the strength of the coupling to the leads and on 
the fraction of potential drop at the contacts. In other 
words, the position of the IETS feature is independent of 
the nature of the contacts. This is the important result 
of this paper, as it means that the IETS signal is more 
stable, in terms of spectroscopic information, than the 
conductance itself. The IETS signal from the internal 
vibrational modes of the molecule is not strongly depen- 
dent on the way the molecule is connected to the leads. 

This occurs because, in spectroscopic terms, the IETS 
signal depends only on the difference between the left 
and right Fermi levels. Thus a vibration mode can only 
be excited by inelastic collision with the charge carriers 
once the bias exceeds the vibrational energy, i.e. V > luq, 
which leads to a corresponding feature in the IETS. Note 
that this is not dependent on the model of the electron- 
vibron coupling — it is purely an effect of non-equilibrium 
inelastic transport at the threshold bias V > Wexti where 
w cx t could be any excitation of the system coupled to the 
injected charge carriers. There are also extreme asym- 
metric cases for which the lineshape of IETS may vary 
or become a dip rather than a peak^, but the position 
of the feature will nonetheless remain at V ~ ujq. 

The conductance peaks at V ~ Sq, however, corre- 
spond to resonant transmission through electronic levels, 
with or without vibration replica^. Such resonances in 
the conductance depend on the relative position of the 
two Fermi levels with respect to the renormalized elec- 
tronic levels, i.e. similar to band offsets in semiconduc- 
tors. The conductance peak position is therefore depen- 
dent on the fraction of potential drop at the contacts. 
In our model, the left Fermi level moves up, while the 
right Fermi level moves down, for positive bias. The 
fraction of how much [i^ moves up for a given bias is 
determined by the factor ijv- The smaller r/y, the larger 
the bias V has to be for hl to become resonant with an 
electronic level sq. Such a mechanism explains the vari- 
ation of the conductance peak positions for the different 
cases (ITV) of fluctuations in the junction that we have 
considered in Fig. [TJ These results confirm experimental 
common knowledge: although an individual molecular 
device based on the use of the conductance as the key 
signal may work wellSS, unless the fabrication of the de- 
vice can be reliably controlled and reproduced then mass 
production will remain impossible. This is especially true 
because the exact nature of the lead-molecule contacts is 
unknown and as yet impossible to reproduce to specifi- 
cation. 

The results we have to chosen to present in this section 
are obtained for a variation of 10% of the hopping inte- 
grals. Obviously, larger variations will lead to stronger 
effects, i.e. more important modification of intensities of 
both the conductance peaks and the IETS features, and 
more important shifts of the position of the conductance 
peak. The IETS features will remain fixed at the same 
bias. It appears here that even small variations of the 



hopping integrals lead to substantial effects. As an il- 
lustration, let us recall that in reality the hopping inte- 
grals between two different electronic orbitals vary expo- 
nentially with the distance between the two atoms sup- 
porting these orbitals. The zero-point motion associated 
with quantum fluctuations are of the order of 0.01-0.05 
A, leading to variations in the hopping integrals of a few 
percent or less. While coordination-induced variations of 
order of 0.1 A in the inter-atomic distance (roughly 10% 
for a carbon double bond, well beyond thermal fluctua- 
tions at room temperature) would lead to changes of 20 
to 50% (and more) of the hopping integrals. Our choice 
of 10% variation of the hopping integrals is intermediate 
between these two regime. 

Furthermore, the analysis performed in this section is 
well suited for large molecular system which have inter- 
nal vibration modes of atoms not involved in the bonding 
of the molecule to the leads. Indeed, for small molecules, 
most, if not all, of the atoms constituting the molecule 
are in close contact with the surface of one or both elec- 
trodes. In these cases the molecular vibronic properties 
(frequency of vibration and electron-vibration coupling 
matrix elements) are strongly dependent on the coupling 
of the molecule to the lead. All (or most) of the atoms in- 
volved in the chemical bonds inside the molecule are also 
involved in the bonding to the leads. These effects have 
been clearly shown in experimental as well as theoretical 
worksM^. 

However for larger molecules (molecular wires) it is 
clear that many vibration modes - the ones mostly lo- 
cated inside the molecular wire - will be much less de- 
pendent on the chemistry of bonding of the end atoms of 
the molecule to the electrodes. These modes are usually 
the optical-type modes and are strongly coupled to the 
LUMO and/or HOMO frontier orbitals as we show in the 
next section. 

It is with these systems in mind that we aim to build 
functionality. Hence the IETS signal of such systenes 
rids us of the need to control accurately the nature of 
the contacts, and is thus a much more useful signal to 
consider when designing and building functional single- 
molecule devices. In order to exploit this phenomenon 
and build a useful device, we need some form of external 
control over the position and/or amplitude of the IETS 
feature. This external control could take several forms 
(magnetic field, chemical concentration, pressure...), but 
here we propose the use of an external electric field, in 
a similar way as a gate voltage is used in Ref. [19| to 
control the conductance. 



III. MODIFYING THE ELECTRONIC AND VIBRONIC 
PROPERTIES OF A MOLECULE WITH A GATE 
POTENTIAL 

We now explore how the properties of the IETS signal 
we have demonstrated in the previous section can be used 
for a more realistic system. 
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For this we use ab-initio calculations to study the ef- 
fect of an uniform electric field, acting as a gate voltage, 
on the electronic and vibronic properties of a realistic 
molecular system. We calculate how the corresponding 
values of loq and 70 in our NEGF calculations are mod- 
ified by the external field for selected vibrational modes 
coupled to the molecular HOMO and lumo levels. We 
then demonstrate in the next section that these results 
can be used to design a selectively functionalized single- 
molecule device. 

As full ab-initio calculations for a realistic single- 
molecule device with three terminals are extremely com- 
putationally demanding, especially for self-consistent cal- 
culations, and are beyond present computational power if 
one has to take the full non-equilibrium and many-body 
effects into account, one has to introduce some approxi- 
mations. 

In the following, as a first step of calculations to 
analyze the potential of functional devices using the 
IETS signal, we perform the calculation for an isolated 
molecule in the presence of an electric field perpendicular 
to the backbone of the molecule. However in real devices 
with metallic electrodes, the electrostatic potential act- 
ing on the molecule from an applied field perpendicular to 
the interelectrode spacing will be substantially distorted, 
on the scale of the molecule, from an uniform field we 
use in the calculations. Hence the results we show in 
this section provide only the general trends of the field 
dependence on the IETS signal for our model ab-initio 
calculations. 

Furthermore, it is known that in real systems the pres- 
ence of the electrode and the coupling of the molecule 
to the electrodes are crucial to determine accurately the 
charge transfer and transport properties of such a molec- 
ular nanojunctions. We have shown in the previous 
section that the nature of the contacts (especially the 
strength of the couplings and the corresponding potential 
drops) do indeed dominate the properties of the conduc- 
tance. The contacts also play an important role in other 
physical transport properties such as heat transport 68 . 

However, we have also shown that the IETS signal is 
virtually independent of these characteristics of the con- 
tacts (hence the study of an isolated molecule). Further- 
more, we concentrate our calculations on the effects of 
the electric field on the vibron modes which have a weak 
amplitudes at the ends of the molecules. Normally, such 
modes would not be strongly coupled to the electrodes 
if the molecule were to be fully connected in a realistic 
molecular device. 

As our test system, we choose the molecule 2,5- 
di [2 '- (para-acetylmercaptophenyl)ethinyl] -4-nitro- 
acetylanilin, shown in Fig. E][a), which has previously 
been used for transport measurements in a break 
junction 3 . Molecules with a similar ethynylphenyl-based 
backbone but with different redox centers in the middle 
benzene ring have also been used in self-assembled 
monolayer transport measurements 6 ^ -. Another reason 
to chose this molecule is based on the fact that the IETS 



signal of conjugated molecules can be modified by a gate 
voltage has been shown in Ref. lfj . Furthermore our 
candidate molecule has peripheral chemical groups that 
provide additional properties. 

This molecule has specific properties that we are able 
to exploit. Firstly, the form of the side-chains provide a 
permanent electric dipole, so that the component of the 
electric field (associated with the gate voltage) that is 
perpendicular to the backbone of the molecule will show 
substantial effects. In particular, the electric field will po- 
larize the electron cloud along the dipole. In the regime 
of strong field, the applied electric field might even bend 
the molecular backbone. Secondly, we will show that the 
dominant phonon modes coupling to the HOMO/lumo 
levels are situated on the central part of the molecule, 
and are thus effectively separated from the leads. 

Because of this, and the properties of the IETS sig- 
nal versus the nature of the contacts, we do not include 
leads in our ab-initio calculations. Rather, we replace 
the end groups of the molecule by hydrogen atoms, and 
then fix these terminal atoms within our supercell when 
structural relaxations are performed. 

We calculate the ground-state electronic and vibronic 
properties of the system using the ABINIT package 60 . 
The calculations are performed with Trouiller-Martin 
pseudopotentials, using the local-density approximation 
(LDA) and the exchange-correlation functional from S. 
Goedecker, M. Teter, and J. Huetter— . We use a super- 
cell of size 50 x 30 x 50 [bohr 3 ] with a plane- wave cutoff of 
30 Ryd at the r-point. The geometries of the molecule in 
the absence and in the presence of the electrical field are 
fully relaxed until the maximum force on each atom is 
less than 0.04 eV/A. These are sufficient for our analysis 
of the frequency variation induced by the external field. 

We first calculate the molecular structure without any 
applied field. The relaxed geometry is planar in the 
(x, y) plane. The molecule has a DFT HOMO-LUMO gap 
of 2.29 eV and presents a permanent dipole moment 
d = do(x, y, 0), lying in the plane of the aromatic cycles, 

of magnitude do = 1-878 D and direction (2, y) — (5, -^) 
where the x-axis is long the molecular backbone. 

The vibronic properties of the molecule are calcu- 
lated from the the dynamical matrix calculated in linear- 
response DFT (as in Ref. [zl): 



D lJ - 



d 2 E t , 



x M,.\l, (h-,.,i)r l , ' 



(3) 



where ri a is the displacement of atom i (of mass Mi) 
a is the reduced direction and E to t is the total energy 
obtained from the DFT calculation. 

The eigenvalues of the dynamical matrix D^b S^ ve ^ ne 
square of the vibron frequencies uj\ while the eigenvec- 
tors are the eigenmodes of vibration V^ a . The electron- 
vibration coupling matrix elements are calculated as 



7fcfc' 




(4) 
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FIG. 4. Molecular wire: (a) ball and stick representation, (b) 
HOMO state, (c) LUMO state, (d — h) Eigenmodes of vibration, 
(d) Mode strongly coupled to the HOMO, mode A = 111 with 
oja ~ 197 meV; (e-f) Modes strongly coupled to the LUMO, 
mode (e) A = 114 with ui\ ~ 206 meV, and mode (f) A = 115 
with lu\ ~ 207 meV; (g-h) Modes strongly coupled to both 
the HOMO and LUMO, mode (g) A = 118 with u)\ ~ 281 meV, 
and mode (h) A = 119 with u\ ~ 284 meV. 
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FIG. 5. Electron-vibration coupling: diagonal matrix ele- 
ments 7fc fc versus the vibration energy uj\ for the two HOMO 
and LUMO Kohn-Sham states k. The amplitude of the matrix 
elements are normalized to arbitrary units, i.e. max(7fc fc ) 2 = 
1. Note that since the IETS signal amplitude around V = ui\ 
is proportional to (7 A ) 2 , the graphs mimics the IETS signal in 
the low-bias region. One would get only two significant peaks 
in the range of applied bias V — to V — 0.35 [V]. 



where H is Hamiltonian of the molecule with atomic po- 
sitions Ti a and \4>k) are the corresponding eigenstates. 

A ball and stick representation of the relaxed molecule, 
the corresponding homo and LUMO states, and the most 
relevant vibrational modes are shown in Fig. |U The cor- 
responding electron-vibration coupling matrix elements 
7fcfc f° r t ne k =homo, LUMO states are shown in Fig. [5] 
Clearly only a few of the vibron modes couple strongly 
to the homo or LUMO states. 

It should also be noted that, at the lowest-order in 
the electron-vibron coupling, the amplitude of the IETS 
features is proportional to (7^) 2 as shown in Refs. [ 
l47l Ir32l . [73|-[76j]. Therefore the graphs in fig. [5] mimic the 
IETS spectrum in the range of applied bias V = to 
V = 0.35 [V]. In that range, one would only see two 
main peaks around V = 0.20 and V = 0.28 [V]. Such 
a behavior justifies a posteriori our NEGF analysis in 
terms of single-mode excitation. 

A central question is how the vibronic properties of 
the molecule (u>\ for selected modes and 7^ fe for the same 
modes and for k = HOMO,LUMO) can be tuned. In our 
case the vibrations will be modified by applying an ex- 
ternal field E to the junction. The external field E acts 
as a potential gate which may control the properties of 
the current flow through the molecule^. 

Calculations of the electronic ground state and of 
the vibronic properties of the molecule are performed 
in the presence of the external electric field using the 
ABINIT package^. Finite electric-field calculations, 
within periodic boundary conditions, are performed by 
introducing an appropriate extra Berry phase in the 
wavefunction£ 7 £~— . 
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FIG. 6. Relative evolution of the vibron energy ui\ and cou- 
pling matrix elements 7 £ fe to the HOMO state as a function of 
the external electric field, for modes A = 111, A = 118 and 
A = 119. The zero-field values for the vibron energies and 
coupling constants (in meV) are uj\ = 197.39 and 7^, =54. 14 
(A = 111); w A =280.74 and 7 £ fe =39.32 (A = 118); oj a =283.93 
and 7^=28.23 (A = 119), respectively. 
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In the present calculations, we take the electric field 
to be lying in the plane of the aromatic cycles and to 
be perpendicular to the main axis x of the molecule, i.e. 
E = Eo(0, 1, 0). In the full three-terminal device, E rep- 
resents the projection of the three-dimensional electric 
field with the plane of the molecule and perpendicular 
to its backbone. This is the component that is expected 
to have the most efficient interaction with the dipole of 
the molecule. We chose that the magnitude of the field 
Eq ranges from 0.001 to 0.020 a.u., i.e. from 0.051 up 
to 1.025 V/A. This values correspond to typical values 
of the electric field between a tip and a sample within a 
scanning probe microscope setup^S. 

Fig. [6] shows how the vibron energy and matrix cou- 
pling element vary, for selected vibron modes, as the ex- 
ternal electric field is increased. We consider the vibron 
modes which are the most strongly coupled to the frontier 
orbitals (Fig. [5). We see that the coupling constant 7^. 
decreases monotonically for modes A=lll and 119 over 
a range of AE = 0.3 [V/A] when E > 0.1 [V/A], with 
an overall decrease of 70% for mode 111, and 80% for 
mode 119. Meanwhile mode 118 decreases much faster, 
showing a switching-like behavior over a smaller range 
AEq — 0.15 [V/A] with an overall reduction of 7^, larger 
than 90%. We also see that the applied external field lifts 
the degeneracy between modes 118 and 119. 

In Fig. [2l we represent how the HOMO and molecular 
level varies with the applied electric field. We see that 
for the values of the field we use, the variations of the 
homo and HOMO-1 levels are only of a few percent. As- 
suming that in a three-terminal device, the Fermi level at 
equilibrium is pined at the mid-gap of the molecule, the 



FIG. 7. Relative evolution of the HOMO and HOMO-1 elec- 
tronic levels e n versus the applied external field. The zero 
field value of the HOMO level is ehomo =-0.18468 Ha = - 
5.0254 eV and e HO MO-i=-0.19478 Ha = -5.3002 eV. The 
top panel shows the relative energy separation between the 
HOMO and HOMO-1 levels Ae = (ehomo — Ew-i)/e with 
£ = \shomo + £homo-i |/2. Two linear regimes or a linear and 
weak quadratic regime are obtained for the dependence of 
ehomo on the external field. The HOMO level varies only by 
a few percent. The applied electric field would not induce 
a transition from the off-resonant to the resonant transport 
regime. Furthermore the energy separation with the HOMO 
and HOMO-1 levels increases with the applied field, validating 
even more the single-level analysis used in our model calcula- 
tions. 



gate voltage is modifying the position of the homo level 
but does not induce a transition between the off-resonant 
transport regime to the resonant transport regime, i.e. 
the homo level is always well below the Fermi level. 
Futhermore 



IV. APPLICATION TO FUNCTIONALITY IN 
SINGLE-MOLECULE DEVICES 

We now concentrate on the feasibility of obtaining a 
functional single-molecule device by using the IETS sig- 
nal. The results of the DFT ab-initio calculations given 
in the previous section will be used as input parameters 
in our NEGF model. The vibron frequency is ujq <— lu\, 
and the electron- vibron coupling constant 7q {"fkk) 2 - 
We take the molecular electronic level Eq to be mid- 
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FIG. 8. Top panel: Evolution of the normalized IETS sig- 
nal around oja = 197 meV for mode A = 111 for the different 
values of external electric field Eo [V/A]. The IETS signal is 
normalized by the value of the current I(V) taken at the same 
bias. The inelastic peak is centered on lo\ = 197 meV for all 
electric field values, only the peak amplitude and background 
changes. Bottom panel: Relative evolution of the normalized 
IETS signal at the inelastic resonance V = uj\ = 197 meV ver- 
sus the values of external electric field Eo- The evolution of 
the signal versus the gate voltage (related to Eo) is typical of 
a non-linear amplification process when Eo decreases from a 
finite value to zero. 

gap, i.e. ±(e LUM o - £homo)A which is ~ ±1.15 eV. 
Finally, we introduce an effective broadening £ol,.r of the 
molecular levels corresponding to the coupling with the 
leads. The coupling is chosen such that the molecular 
level broadening is less than the spacing between levels 
(ehomo — ehomo-i )/3.3 ~ 80 meV. 

Using these values, we calculate the IETS signal 
around V ~ u)o, and study how this signal is modified 
by the external applied field E = Eq(0, 1,0). The up- 
per panels of Figs. \E\ and O show the evolution of the the 
IETS signal for NEGF calculations with the parameters 
for mode 111 and mode 118. Here we have chosen to nor- 
malize the IETS by the current itself, as it reduces the 
effect of the slope background on the IETS peaks (see 
Fig.©. 

The lower panels of Figs. [5] and [9] show the relative evo- 



FIG. 9. Top panel: Evolution of the normalized IETS signal 
around u;a=280 meV for mode A = 118 for the different values 
of external electric field Eo- Bottom panel: Relative evolution 
of the normalized IETS signal at V = cj\=280 meV versus the 
values of external electric field Eo. The evolution of the signal 
versus the gate voltage (related to Eo) is typical of an strong 
switching process from high to low values when Eo increases 
from 0.05 to 0.20 [V/A]. 



lution of the normalized IETS signal for a fixed source- 
drain applied bias versus the external electric field act- 
ing as a gate voltage. These results demonstrate that we 
can indeed achieve a functional behavior from the IETS 
signal when working at low applied bias around the fre- 
quencies of the vibron modes (here modes 111 and 118) 
that are strongly coupled to the HOMO molecular level. 

The behavior of the IETS signal in Fig. [5] is typical 
of a non-linear amplification of the source-drain current 
versus the gate voltage represented by Eg. The signal 
is amplified in a non-linear manner when Eq decreases 
from a finite value to zero. Note that the curve is simply 
convex with no inflection point. 

In Fig. [9] we obtain a typical switching behavior of the 
IETS signal from two plateaux of high and low value as 
Eq increases. Note that now the curve presents an inflec- 
tion point around i?o=0-10 [V/A]; and that the switch- 
ing occurs over a small range of applied field Ai?o=0.15 

[v/A]. 
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The relative variations of the current-normalized IETS 
signal is of the order of 20-25% and is not as important 
as the relative variation of the corresponding coupling 
matrix elements (j^ k ) 2 shown in Fig. [5] This is to be ex- 
pected, since the full non-equilibrium inelastic transport 
properties are not simply related in a linear way to the 
square of the coupling matrix elements (t^) 2 , as would 
be obtained from perturbation theory 81 . In a full non- 
equilibrium calculation, complex non-linear effects enter 
into accoun t 40 ' 57 and cannot be described perturbatively. 

In terms of practical devices, the IETS signal could 
be measured by using an electronic circuit similar to 
that developed for non-linear amplification through sin- 
gle fullerene molecules 8 ^. However, in order to extract 
the IETS, one would use a lock-in technique for the in- 
put source-drain bias, and measure the second harmonic 
of output voltage at the load resistance 9 -. One then gets 
a signal proportional to the IETS signal since the device 
works within the linear regime for the source-drain cur- 
rent versus source-drain bias for low applied bias around 
the vibron frequencies (100-350 meV). 

Until now, IETS measurements in molecu- 
lar nanojunctions have been performed at low 
temperature 4 -' 1 ^ 14 ' 16- — . This is necessary in order 
to have good mechanical stability of the molecular 
junctions, as well as a good resolution of the IETS 
feature which depends on the thermal fluctuations of 
the soft acoustic-like vibration modes 8 -. The width of 
the IETS features increases with increasing tempera- 
ture, with a loss of resolution above liquid nitrogen 
temperatur e) 12 ' 19 . All this limits the use of the single- 
molecule devices to temperatures below ~70K. It can 
be expected that by working with other candidates of 
molecular connectors, or by using different concepts to 
build three-terminal devices (i.e. for example by using 
carbon-based electrodes 8 ^), this temperature limit could 
be lifted. 



V. DISCUSSION 

Using a two-step theoretical framework we have shown 
that the IETS signal through molecular junctions can 
be used to achieve functionality within a single-molecule 
device. By using a NEGF approach for a model sys- 
tem, we have shown that the IETS is virtually indepen- 
dent of the nature of the contact between the molecule 
and the source and drain. The electronic and vibronic 
properties of a realistic molecular candidate (with an 
ethynylphenyl-based backbone) are calculated using an 
ab-initio method. We have shown how an externally ap- 
plied electric field, simulating the presence of a gate elec- 
trode, can control the vibron properties of the molecule 
and therefore the corresponding IETS signal. 

Multi-functionality is demonstrated within a single 
molecule: non-linear amplification and switching are 
both present for different vibron modes in the molecu- 
lar system we have studied. Such functionality should 



be reproducible from device to device, since the IETS is 
virtually independent of variations in the molecule-lead 
contacts. 

Furthermore, a recent theoretical study on the diffi- 
culty of gate control in molecular transistors 8 - 5 - has shown 
that when the molecular energy levels are away from the 
Fermi level (i.e. the off-resonant transport regime), they 
can be shifted by the gate voltage. However when the 
molecular levels move close to the Fermi level, the shifts 
become extremely small and almost independent of the 
gate voltage. This indicates that it may be difficult to use 
the gate voltage to control transition in the conductance 
between the off- resonant regime ("OFF" state) and the 
resonant regime ("ON" state). This difficulty does not 
occur in our scheme, since we do not rely on the use of 
the conductance as the functional signal. 

In the present study we have considered the low 
electric-field regime for the control of the IETS at finite 
bias. In this regime, most of the functionality comes from 
the dependence of the coupling matrix elements 7^. on 
the electric field, while the vibron frequencies are more 
or less constant. In the regime of stronger electric field, 
the polarization of the electronic clouds may be strong 
enough to lead to modification of some chemical bonds in 
the molecule. We expect then that the variation of the vi- 
bron frequencies will be important and the displacement 
of the vibron frequencies will dominate the functionality 
of the single-molecule device. 

We can now suggest designing single-molecule sensors 
in the following manner: one uses a molecule with pe- 
ripheral chemical groups that can actively react with sur- 
rounding molecules. After chemical reaction, the vibron 
modes of the molecular backbone that couple strongly 
to the HOMO-LUMO molecular states, and/or the corre- 
sponding coupling matrix elements, are modified. One 
can then monitor the state of the molecular sensor by 
measuring the evolution of the IETS signal. Finally, the 
mechanisms we describe above are most relevant in the 
vibrational theory of the sense of smell (olfaction)^— . 
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